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We present a new measurement of air shower muons made 
during atmospheric ascent of the High Energy Antimatter 
Telescope balloon experiment. The muon charge ratio jj,'^ 
is presented as a function of atmospheric depth in the momen- 
tum interval 0.3-0.9 GeV/c. The differential momentum 
spectra are presented between 0.3 and ~50 GeV/c at atmo- 
spheric depths between 13 and 960 g/cm^. We compare our 
measurements with other recent data and with Monte Carlo 
calculations of the same type as those used in predicting at- 
mospheric neutrino fluxes. We find that our measured fi" 
fluxes are smaller than the predictions by as much as 70% at 
shallow atmospheric depths, by ~20% at the depth of shower 
maximum, and are in good agreement with the predictions at 
greater depths. We explore the consequences of this on the 
question of atmospheric neutrino production. 

PACS numbers: 96.40. Tv, 14.60.Pq, 14.60.Ef 



I. INTRODUCTION 

Measurements of the flux of atmospheric neutrinos by 
large underground detectors have consistently disagreed 
with theoretical predictions, a fact that has been inter- 
preted in terms of possible neutrino oscillations. The 
most compelling evidence thus far comes from the Su- 
perKamiokande experiment The discrepancies be- 
tween experiment and theory are well beyond the statisti- 
cal uncertainties in the measurements. The correct inter- 
pretation of the effect requires a detailed understanding 
of neutrino production in the atmosphere. 



Model predictions of the absolute flux of neutrinos 
produced in air showers are uncertain due to a number 
of difficulties: one problem is the normalization of the 
primary cosmic-ray flux, for which measurements vary 
by ±15% or more. Because of this, different assumptions 
are made by different authors. For instance, the model 
of Honda et al. Q , used as the starting point in the anal- 
yses of the SuperKamiokande collaboration, assumes a 
primary cosmic ray flux normalization that is in excess 
of the current measurements by balloon experiments [Q, 
by as much as ~20% in the energy range beyond about 
5 GeV, which is relevant to the production of sub- to 
multi-GeV atmospheric neutrinos. In contrast, the model 
of the Bartol group uses a primary spectrum in good 
agreement with the most recent measurements. Different 
assumptions are also made about the details of particle 
production in atmospheric interactions, and about the 
treatment of geomagnetic effects (see Ref. ^). Thus, 
while the authors of these calculations expect that the 
predictions of absolute neutrino fluxes might have a to- 
tal uncertainty of about 18% it seems fortuitous that 
the predicted neutrino rates have been found to agree 
with each other within about 10%. It should also be 
noted that in the neutrino oscillation analysis of the Su- 
perKamiokande collaboration, an additional normaliza- 
tion factor of 1.16 in the cosmic ray intensity is intro- 
duced to obtain the best oscillation fit to the measured 
muon and electron neutrino rates as a function of zenith 
angle. 

Neutrino production in the atmosphere is closely cou- 
pled to muon production, as both types of particles are 
produced together in pion and kaon decays, and as some 
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of the muons decay to produce further neutrinos. Monte 
Carlo simulations of neutrino production naturally pre- 
dict the spectra of muons as a function of atmospheric 
depth. Therefore, a detailed experimental test of the 
predictions is possible through measurements of the in- 
tensity of muons at different depths in the atmosphere, 
for instance with a balloon-borne particle detector during 
atmospheric ascent. If such a detector includes a mag- 
net spectrometer, separate measurements of the fj,~ and 
/i"*" fluxes are possible. This is relatively easy for neg- 
ative muons, as negatively charged particles other than 
electrons are rare in air showers (and electrons are easily 
rejected), while non-interacting protons generate a large 
background for positive muons above ~1 GeV. Such mea- 
surements have been made by a number of instruments, 
including the MASS |lO0, HEAT [^,[3), CAPRICE 
and IMAX |jl^ magnet spectrometers and older, 
less sensitive balloon payloads. 

Here we describe measurements made with the High 
Energy Antimatter Telescope (HEAT), a balloon instru- 
ment optimized for the study of high-energy cosmic-ray 
electrons and positrons. HEAT uses several complemen- 
tary particle identification techniques, which are also well 
suited to the identification of muons during atmospheric 
ascent of the balloon. This instrument was flown twice, 
in 1994 and 1995. Preliminary results from the first 
flight on relative abundances of /i^ and /i^ were reported 
previously [p^ . The present paper describes measure- 
ments made during the second balloon flight: the rela- 
tive abundances of fj,~^ and fi~ at momenta between 0.3 
and 0.9 GeV/c, and the differential spectra of /i~ at mo- 
menta between 0.3 and 50 GeV/c, at atmospheric depths 
between 13 and 960 g/cm^. We compare the results with 
other measurements and with calculations made with a 
modified version of the TARGET Monte Carlo algorithm 

, developed by the Bartol group to predict atmospheric 
neutrino rates. 



II. MUON IDENTIFICATION AND 
BACKGROUNDS 

The HEAT instrument is described in detail elsewhere 
p6| . It combines a superconducting magnet spectrom- 
eter (using a drift-tube tracking hodoscope), a time- 
of-flight (TOF) system, a transition-radiation detector 
(TRD) and an electromagnetic calorimeter (EMC). For 
its second flight, it was launched from Lynn Lake, Man- 
itoba, Canada, at a vertical geomagnetic cutoff rigidity 
of well below 1 GV. The flight took place on 23 Au- 
gust 1995, approximately at minimum solar activity. The 
ascent from an atmospheric overburden at the ground 
of 960 g/cm^ to a float altitude at 3 g/cm^ required 
3.1 hours, during which charged atmospheric secondary 
particles were detected and recorded. Near the end of 
ascent, at about 13 g/cm^, the instrument trigger config- 
uration was changed to commence preferential measure- 
ments of electrons and positrons, affecting the ability to 



measure absolute muon intensities at shallower depths. 
Therefore we report measurements of the muon charge 
ratio at all atmospheric depths, but measurements of ab- 
solute muon intensities only at depths greater than 13 
g/cm^. 

In identifying muon events, TRD information is not 
used. The TOF system measures the velocity of the 
particle /3 = v/c with a resolution of ct/j = 0.15, per- 
mitting complete rejection of upward-going particles; in 
addition, the amount of light generated in the scintilla- 
tion counters of the TOF system is used to measure the 
magnitude of the particle's electric charge Ze with a res- 
olution of az = 0.11, permitting the unambiguous iden- 
tification of singly-charged particles. The magnet spec- 
trometer measures the sign of the particle's charge from 
the direction of the deflection in a magnetic field of about 
1 T, with a field integral over the particle's trajectory of 
J B ■ dl ^ 4.2 kGm. The magnetic rigidity R = pc/Ze 
of the particle is determined from the amount of deflec- 
tion; at rigidities between 0.3 and 0.9 GV, the resolution 
achieved is an = 0.08 to 0.11 GV. The EMC records the 
pattern of energy deposits in 10 scintillation counters, 
each topped by a 0.9 radiation length-thick lead sheet. 
In each layer, the energy deposited is measured in units 
of the energy loss by a vertical minimum-ionizing particle 
(MIP). A shower sum is obtained by adding the signals 
from the 10 scintillators. 

The selection of muon events proceeds in three steps. 
First, a high-quality spectrometer track and measure- 
ment of the rigidity R are required, as described in 
Ref. ||l^. Also required are down-going and singly- 
charged TOF characteristics. Second, as shown in Fig- 
ure |l], the distribution of reconstructed particle rigidity 
R is studied as a function of velocity /3, and compared 
with ideal curves R = (to/Z)/3/V(1 - P^) where m is the 
mass and Z the charge of the particle. Clearly identifiable 
in the figure are populations due to He (and d), p, /i^ 
and iJ,~ events, tt* and some events cannot be distin- 
guished from the /i* populations, and are a small back- 
ground to the muon signal. We select /i-like events by 
requiring: 1) (3 > 0.85 and 2) 0.3 < \R\ < 0.9 GV, for the 
low-energy ratio, and 2) R < —0.3 GV for the ii~ 

energy spectra. For R > 0.9 GV, fj,~^ events become in- 
distinguishable from non-showering hadrons (mostly pro- 
tons), so that /z"*" spectra are not measured with this 
instrument. Third, we study the behavior of the shower 
sum measured by the EMC as a function of particle rigid- 
ity R, as shown in Figure |^. Ideal curves for electrons 
and positrons are obtained by assuming a simple linear 
relationship between EMC sum (proportional to energy) 
and R. For heavier particles, a calculation of EMC sum 
is made by integrating Bethe-Bloch energy losses within 
the EMC scintillators. Populations of events due to e^, 
p and /i^/tt^/A'^ are identifiable. By augmenting the 
selection criteria of Figure |^ with the requirement that 
particles not shower in the EMC (EMC sum < 15), e='= 
events with \R\ > 0.3 GV are rejected. The low-rigidity 
proton events that range out in the EMC and would ap- 
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pear to contaminate the /z"*" population are rejected by 
the /3 > 0.85 requirement of Figure |l|. 

With these criteria, we achieve essentially complete 
rejection of electron events, but there remains a small 
background to the muon signal due to pions and kaons. 
We estimate from Monte Carlo simulations of air show- 
ers based on CERN's GEANT-FLUKA algorithms 
that the tt^ flux at a depth of 13 g/cm^ is only 2% that 
of , in agreement with another calculation jlSj - with 
fluxes at a much lower level - and that this fur- 
ther decreases with increasing atmospheric depth. More- 
over, only pions that do not interact can be mistaken for 
muons, which occurs 39% of the time, so that the back- 
ground to the muon measurement due to atmospheric 
pions is only about 0.7-0.9% near float, decreasing to less 
than 0.4% at depths greater than 300 g/cm^. Such a 
small background is not corrected for here. Occasionally, 
cosmic-ray interactions in the instrument result in tt^ 
production, at an even more modest level than the at- 
mospheric pion background. GEANT-based simulations 
indicate that only 0.04% of proton-induced events yield 
a misidentification as a As the pay load slowly ro- 
tated throughout the flight (as determined with a solar 
sensor attached outside the gondola), and did not align 
itself with the Earth's magnetic field, possible geomag- 
netic East- West asymmetries are averaged out. Further- 
more, such asymmetries in the primary proton flux are 
only expected at momenta near the geomagnetic cutoff, 
which for the Lynn Lake flight is well below the energies 
of interest here. 



III. RESULTS AND DISCUSSION 

A. Muon Charge Ratio 

The number of low-energy fj,~^ and fi~ events detected 
as a function of atmospheric depth and the ^'^ / fj,~ ratio 
are summarized in Table || The jX^ / ^~ ratio is shown 
as a function of atmospheric depth in Figure || (labeled 
HEAT 95), together with the measured ratio from the 
first HEAT flight |l|l_(labeled HEAT 94) and other re- 
cent measurements ||14p^ , H - As noted on the figure, 
the various analyses have used different magnetic rigid- 
ity ranges, and moreover, the measurements were made 
at different solar epochs and different geomagnetic rigid- 
ity cutoffs Rcutoff , so that the data are not truly directly 
comparable. The two HEAT measurements of the muon 
ratio are the most statistically significant, and are es- 
sentially consistent with each other and with other mea- 
surements, within the appreciable errors. Based on the 
HEAT measurements, no clear correlation of the muon 
ratio with geomagnetic cutoff rigidity is observed. The 
depth dependence of the charge ratio seems to be essen- 
tially flat, although a slight decrease from about 1.3-1.4 
at high altitudes (3-50 g/cm^) down to about 1.1 at the 
ground cannot be excluded. We note that the charge 



ratio measured by both HEAT 95 and CAPRICE 94 at 
small depths appears anomalously high. This effect is not 
understood at present. Also shown on Figure^ are calcu- 
lations with the TARGET algorithm |^ (widely used for 
neutrino flux calculations), for conditions of solar min- 
imum and maximum activity. Both HEAT flights oc- 
curred under essentially solar-minimum conditions. The 
TARGET calculations have been made for average pri- 
mary fluxes at solar minimum and maximum, which may 
not exactly represent the actual spectrum at the time 
of the flight. The calculations are for a location with 
no geomagnetic cutoff, and are intended for comparison 
with the HEAT 95 data only. The agreement between 
the HEAT measurements and the solar-minimum TAR- 
GET calculation appears to be fairly good. Note that 
the fluctuations in the simulated distributions at shallow 
atmospheric depths are statistical, and indicative of the 
small number of muons having been produced in the air 
showers at the highest altitudes. 

B. Energy Spectra of Negative Muons 

The absolute intensity of in a rigidity interval Ai? 
with an average rigidity i?, at atmospheric depth d is 
obtained with: 



^scan^acc 

where N^^ is the number of events recorded at (d, R), 
At is the time spent at depth d, is the live time frac- 
tion, 6/ is an event-transmission loss correction, eg is a 
correction to account for the fact that the muon flux is 
increasingly less isotropic deeper in the atmosphere, e^t 
is the efficiency of the basic event cleanliness criteria ap- 
plied to the drift-tube hodoscope track, (JIA) is the ge- 
ometrical factor, escan is a "scanning efficiency" correc- 
tion (described below), and eacc is the muon acceptance 
efficiency. The acceptance of the instrument decreases 
rapidly for particles incident at a zenith angle greater 
than about 25°, so that particle intensities reported here 
are essentially for vertical incidence. 

Both R and AR are weighted to account for the details 
of the energy spectrum, according to: 

_ ^ Rf{R)dR ^ f{R)dR 

!^:f{R)dR f(R) 

where f{R) oc i?"" is the rigidity power-law spectrum, 
with spectral index a varying between —0.56 and 3.5 
depending on both R and d {a is experimentally de- 
termined from the spectra before any of the normaliza- 
tion corrections are applied). At is measured with an 
on-board clock, e* is determined using on-board scalers 
which count clock cycles while the instrument is available 
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for a trigger or busy processing an event, and e; is deter- 
mined by careful accounting of event numbers generated 
on board compared to events successfully transmitted. 
eg is calculated using a standard prescription po|] , where 
the zenith dependence of the muon flux is taken to be 
PQgn(d) -p^fith the exponent n a function of atmospheric 
depth d; we have used n{d) — (d/lOSOg/cm^) x Useaievei 
with Useaievei =2. e^t is obtained by a careful account- 
ing of the number of events recorded compared to the 
number of events with a successful minimal track recon- 
struction. (nA) and Cqcc are determined with the aid of 
a GEANT-based simulation of the response of the HEAT 
instrument. Cscan is a correction factor introduced based 
on the visual scanning of several hundred events to ac- 
count for residual differences between the reconstruction 
efficiency of real events compared to that of simulated 
ones, and is found to be escan = (0.9 ± 0.1). The various 
parameters described above are given in Tables || and 

The final intensities as a function of momentum for 
various atmospheric depths are shown in Figure ^ and 
given in Table IV . Also shown in the figure are the mea- 
surements of the MASS 0,|lltl and CAPRICE |l| exper- 
iments. The HEAT sample of 10327 events collected 
during ascent is to be compared with the MASS samples 
of 2893 events (1989 flight) and 4471 events (1991 flight) 
and the CAPRICE sample of 4627 events. Although the 
1989 MASS measurements were also made in Northern 
Canada (from Prince Albert, Saskatchewan), the flight 
occurred at a different solar epoch (1989), at the time of 
a significant Forbush decrease. The 1991 MASS measure- 
ments were made from Fort Sumner, New Mexico. The 
CAPRICE data were collected at Lynn Lake, in 1994, 
and so are more directly comparable with our measure- 
ments. The general level of agreement between the data 
sets should be noted. 



C. Comparison with Model Calculations 

1. One- Dimensional TARGET Algorithm 

In Figure ^, we compare the HEAT measurements re- 
orted here with predictions of the TARGET algorithm 
for conditions of solar minimum and maximum activ- 
ity, shown as solid and dotted curves, respectively. (The 
solar-minimum curves are the ones of interest here, but 
the solar maximum curves are shown as well to illustrate 
the extent of the effect of the solar cycle on muon pro- 
duction.) These curves are obtained with the standard 
TARGET algorithm, which simulates vertically incident 
cosmic rays, and which follows the development of the 
air shower in one dimension only. No corrections for ge- 
omagnetic effects are made. This is the algorithm devel- 
oped by the Bartol group and used in predictingunder- 
ground neutrino rates from atmospheric sources |||,^. In 
Figure o, we show the /i" growth curves for different mo- 



mentum intervals, also compared with the ID TARGET- 
based predictions for solar minimum conditions (solid 
curves). The calculations were not made for the high- 
est momentum bin. 

There is general similarity between the experimental 
and simulated distributions, with some notable differ- 
ences. For instance, the predictions are significantly in 
excess of the measurements below 4 GeV/c at atmo- 
spheric depths between 13 and 250 g/cvc?. The ratio 
of simulated to measured intensity varies from 1.2 ± 0.2 
near shower maximum at 200 g/cm^ to 1.7±0.3 at depths 
between 13 and 140 g/cm^. At depths beyond shower 
maximum, or at momenta greater than 4 GeV/c, the 
simulations agree very well with our measurements. A 
similar trend was found by the CAPRICE collaboration 
Jl^ : simulations predict more /i~ events than they mea- 
sure below about 1 GeV/c, but they find that the ratio 
of simulated to measured intensity is greatest at shower 
maximum, with a value of 1.8 ± 0.1. 



2. Three-Dimensional TARGET Algorithm 

In an attempt to understand the origin of the discrep- 
ancy between the predicted and measured muon intensi- 
ties, a new version of the TARGET algorithm was pro- 
duced in collaboration with T. Gaisser and T. Stanev 
of the Bartol Research Institute. In this, three dimen- 
sional air shower development effects are taken into ac- 
count, and the primary cosmic ray arrival direction is 
sampled isotropically, rather than assuming vertical in- 
cidence. Geomagnetic effects are not yet included in the 
calculations. Figure ^ also shows the /z~ momentum 
spectra at various atmospheric depths obtained with the 
3-dimensional TARGET algorithm (dashed curves), for 
solar minimum conditions. The 3D calculations are in 
substantially better agreement with the data than the 
ID calculations. In Figure ^, the measured growth 
curves are also compared with 3D TARGET predictions 
(dashed curves). Here again the 3D predictions are a 
more adequate representation of the data. 

Figure |^ shows distributions of the ratios of predicted 
to measured /i~ intensities, for ID TARGET calcula- 
tions (top panel) and 3D TARGET calculations (bottom 
panel), respectively. These are cumulative distributions 
for all momentum and atmospheric depth bins. Each 
ratio is weighted by the square of the error on the ra- 
tio derived from the experimental error on the intensity. 
The distribution for ID calculations has a mean 1.13, in- 
dicating an average overprediction of 13%, whereas for 
3D calculations the mean is 1.07, a slightly better agree- 
ment. The main improvement however is in the reduced 
RMS variance of 0.17 for the distribution for 3D calcula- 
tions compared to 0.27 for the ID calculations. Thus, the 
more realistic calculations that take into account 3D air 
shower development and primary zenith arrival direction 
constitute a clear improvement in the representation of 
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muon production. 

3. Neutrino Production 

The ID and 3D TARGET algorithms were used to pre- 
dict (Vfj, +T'fj^) and (i/o+Fo) intensities at different atmo- 
spheric depths and in different momentum bins. The 
calculations are made for solar minimum conditions, for 
no geomagnetic rigidity cutoff, and for primary cosmic 
rays arriving within 30° of the zenith. The resulting neu- 
trino growth curves in different momentum intervals are 
shown in Figure ^ The calculations are made only up 
to 32 GeV/c. Although there are differences between 
the ID and 3D predictions at momenta less than about 
1 GeV/c, these differences are most important at mid- 
to-high altitudes. The neutrino intensities at the ground 
level, which are the ones of relevance to the underground 
neutrino studies, are summarized in Table [v| for ID and 
3D calculations. The neutrino intensities at ground level 
appear not to be altered much by the 3D effects. A simi- 
lar conclusion was also reached in a study by Battistoni et 
al. |2l] , where detailed calculations of atmospheric muon 
and neutrino production are made in one and three di- 
mensions. 



IV. CONCLUSIONS 

We have made statistically significant measurements of 
air shower muons as a function of atmospheric depth. We 
report the muon charge ratio fj,~^ / ^~ in the momentum 
range 0.3-0.9 GeV/c and the momentum spectra of fi^ 
in the range 0.3-50 GeV/c, at atmospheric depths from 
f3 to 960 g/cm^. The charge ratio is essentially con- 
stant with altitude within errors, with a possible decrease 
from f.3-1.4 at high altitudes to f.l at the ground. A 
comparison of our measured /i~ momentum distributions 
with model calculations indicates significant discrepan- 
cies with the predictions of the standard one-dimensional 
TARGET algorithm: our measured fiuxes are lower than 
the calculated ones at shallow depths before about shower 
maximum. Calculations of the muon intensities with a 
new version of the TARGET algorithm, accounting for 
three-dimensional air shower development, lead to a sub- 
stantially improved agreement with our data. A detailed 
representation of atmospheric secondary production thus 
benefits from the more realistic simulations. The aver- 
age excess of about 7% of the 3D calculations over our 
measured intensities is comparable to possible systematic 
effects in our experiment. Thus, within this uncertainty, 
the 3D TARGET algorithm generates atmospheric sec- 
ondary particle intensities which are in agreement with 
the measurements. 

The three-dimensional air shower development effects 
do not appear to impact significantly the atmospheric 
neutrino rates at the ground, but merely the pattern of 



neutrino production altitudes. Thus, we estimate that 
the neutrino intensities predicted by the ID version of the 
TARGET algorithm are also accurate to about 7%. This 
is to be compared with the accuracy of 14-18% first esti- 
mated for such calculations Even though the Honda 
et al. model uses different assumptions about the pri- 
mary cosmic ray flux and about the atmospheric inter- 
action characteristics, it predicts neutrino intensities on 
the ground which agree with the TARGET predictions 
and with our data at the level of 7-10%. Therefore, one 
might put into question the additional normalization fac- 
tor of 1.16 of the primary cosmic ray spectrum that is 
introduced in the SuperKamiokande neutrino oscillation 
analysis. 
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FIG. 1. Distribution of particle rigidity 7? as a function of 
velocity {3 = v/c for ascent data. 
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FIG. 2. Distribution of shower sum in the EMC as a func- 
tion of particle rigidity R for ascent data. 
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FIG. 3. Measured muon ratio /.i^ / ^ as a function of atmo- 
spheric depth. The curves are calculations with the TARGET 
algorithm. The HEAT data are for 0.3 < |-R| < 0.9 GV. 
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FIG. 4. Differential n~ momentum spectra as a function of 

atmospheric depth. The HEAT measurements axe compared 

with the MASS and CAPRICE measurements for quaUtative 

purposes only. 



TARGET ID solar min TARGET 3D solar min 

TARGET ID solar max 



O -2 

^. 10 



,10 



" -4 

.2 10 




10 



X 



'10 



10 



10 



10 



10 



HEAT 95 13-32 g/cm 



HEAT 95 67-140 g/cm 




• HEAT 95 250-350 g/cm' 




1 10 



10 



• HEAT 95 32-67 g/cm 




HEAT 95 140-250 g/cm 




• HEAT 95 350-840 g/cm 




• HEAT 95 840-960 g/cm' 




1 10 
Momentum (GeV/c) 

FIG. 5. Differential n~ momentum spectra measured by 
HEAT for various atmospheric depths. The curves are predic- 
tions of the ID TARGET algorithm at solar minimum (solid 

curves) or solar maximum (dotted curves) conditions, and the 
3D TARGET algorithm at solar minimum (dashed curves), 
respectively. 
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FIG. 6. n~ growth curves measured by HEAT for different 
momentum intervals. The curves are predictions of the ID 
TARGET algorithm (soUd curves) and 3D TARGET algo- 
rithm (dashed curves), respectively, for solar minimum con- 
ditions. 
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FIG. 7. Distribution of the ratios of predicted to measured 

fi~ intensities, for ID TARGET calculations (top panel) and 
3D TARGET calculations (bottom panel), respectively. Each 
ratio is weighted by the square of the error on the ratio derived 
from the experimental error on the intensity. 
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FIG. 8. u growth curves predicted by TARGET at Lynn 
Lake, for solar minimum conditions, in different momentum 
intervals. The curves are: ID (I'p + i'^), solid curves; 3D 
{v^ + Vfi), dashed curves; ID (z/e + 1^6), dotted curves; 3D 
{ve+Ve), dot-dashed curves. 
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TABLE I. Muon counts at 0.3 < \R\ < 0.9 GV as a func- 
tion of atmospheric depth d. Uncertainties are statistical. 



d 

(g/cm^) 


d 

(g/cm^) 








3-4 


3.45 


134 


65 


2.06 ±0.31 


4-7 


5.13 


35 


27 


1.30 ±0.33 


7-13 


9.85 


53 


37 


1.43 ±0.31 


13-32 


23.2 


305 


193 


1.58 ±0.15 


32-67 


49.1 


458 


348 


1.316 ±0.094 


67-140 


105 


1264 


954 


1.325 ±0.057 


140-250 


190 


1691 


1463 


1.156 ±0.041 


250-350 


298 


605 


540 


1.120 ± 0.066 


350-840 


499 


810 


792 


1.023 ± 0.051 


840-960 


957 


1076 


990 


1.087 ± 0.048 



TABLE IL Rigidity-dependent parameters and corrections 
used in the determination of the absolute /i~ intensities. 



R 


R 




(flA) 


^acc 


(GV) 


(GV) 


(%) 


(cm'^sr) 


(%) 


0.3-0.5 


0.40 


68.35±0.22 


568.9±2.6 


52.0±0.5 


0.5-0.8 


0.65 


74.46±0.20 


613.1±2.8 


63.0±0.5 


0.8-1.2 


0.99 


78.77±0.19 


608.2±2.7 


66.1±0.5 


1.2-1.7 


1.43 


80.46±0.18 


604.4±2.7 


65.4±0.5 


1.7-2.5 


2.06 


81.57±0.18 


600.0±2.7 


64.6±0.5 


2.5-4 


3.13 


81.64±0.18 


601.9±2.7 


63.9±0.5 


4-8 


5.52 


81.48±0.18 


598.9±2.7 


62.5±0.5 


8-16 


11.0 


81.45±0.18 


603.6±2.7 


61.2±0.5 


16-32 


21.3 


81.45±0.18 


606.0±2.7 


47.0±0.4 


32-50 


39.2 


81.43±0.18 


602.4±2.7 


15.5±0.2 



TABLE III. Dopth-dcpondont parameters and corrections 
used in the determination of the absolute jj," intensities. 



d 


d 


At 


£t 


ei 


ee 


(g/cm^) 


(g/cm^) 


(±ls) 


(%) 


(%) 


(%) 


13-32 


23.2 


1453 


48.1±1.8 


89.190±0.079 


99.68 


32-67 


49.1 


1307 


41.9±1.6 


96.585±0.047 


99.33 


67-140 


105 


1571 


49.5±4.6 


95.051±0.051 


98.58 


140-250 


190 


1285 


65.8±6.9 


92.894±0.072 


97.45 


250-350 


298 


484 


82.9±0.3 


95.45±0.12 


96.05 


350-840 


499 


1163 


95.7±0.2 


94.30±0.14 


93.50 


840-960 


957 


5335 


99.3±0.1 


97.09±0.10 


88.03 
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TABLE IV. Number of recorded /U events and their intensity (in (cm^s sr GeV/c) ^) as a function of momentum and 
atmospheric depth. 





0.40 GeV/c 


0.65 GeV/c 


0.99 GeV/c 


1.43 GeV/c 


2.06 GeV/c 


23.2 


81 


86 


89 


42 


51 


g/cm 


(4.12 ± 0.67) X 10" 


"'^ (2.06 ± 0.33) X 10"^ 


(1.45 ± 0.23) X 10"^ 


(5.5 ± 1.1) X 10~* 


(4.22 ± 0.77) X 10"" 


49.1 


135 


181 


100 


97 


90 


g/cm 


(6.7 ± 1.0) X 10"- 


^ (4.42 ± 0.62) X 10"^ 


(1.68 ±0.26) X 10"^ 


(1.29 ±0.20) X 10"^ 


(7.6 ± 1.2) X 10"" 


105 


354 


494 


403 


285 


241 


g/cm^ 


(1.32 ± 0.20) X 10" 


"^ (8.7 ± 1.3) X 10"^ 


(4.86 ± 0.75) X 10"^ 


(2.74 ± 0.43) X 10"^ 


(1.46 ± 0.23) X 10"^ 


190 


568 


721 


626 


441 


370 


g/cm 




-2 Cl 91 + IQ^I V in~^ 


(7 O + T -\\ y -\ 


C4 + n v in~^ 

l^.OO \J.\JOJ A ±\J 


I^.j-O in u.o^j A ±u 


298 


198 


268 


245 


191 


189 


g/cm^ 


(1.47 ±0.20) X 10" 


"^ (9.4 ± 1.3) X 10"^ 


(5.84 ± 0.79) X 10"^ 


(3.63 ±0.51) X 10"^ 


(2.26 ±0.31) X 10"^ 


499 


235 


441 


398 


309 


275 


g/cm^ 


(6.68 ±0.89) X 10" 


(5.83 ± 0.73) X 10""^ 


(3.56 ±0.45) X 10"^ 


(2.21 ±0.29) X 10"^ 


(1.23 ±0.16) X 10"^ 


957 


280 


550 


624 


630 


653 


g/cm^ 


(1.85 ± 0.23) X 10" 


"^ (1.61 ±0.19) X 10"^ 


(1.22 ± 0.14) X 10"^ 


(9.8 ± 1.2) X 10"" 


(6.36 ± 0.75) X 10"" 






3.13 GeV/c 


5.52 GeV/c 


11.0 GeV/c 


21.3 GeV/c 


39.2 GeV/c 


23.2 


34 


23 


7 


4 





g/cm^ 


(1.52 ±0.32) X 10 


(4.04 ± 0.97) X 10"^ 


(6.5 ±2.6) X 10"® 


(2.4 ±1.2) X 10"® 




49.1 


62 


44 


22 


4 


1 


g/cm^ 


(2.80 ± 0.49) X 10 


(8.0 ± 1.5) X 10-^5 


(2.02 ± 0.49) X 10"® 


(2.6 ± 1.3) X 10"® 


(1.8 ± 1.8) X 10"® 


105 


173 


128 


44 


12 


3 


g/cm^ 


(5.73 ± 0.94) X 10 


(1.68 ±0.29) X 10"* 


(2.93 ± 0.61) X 10"® 


(5.3 ± 1.7) X 10"® 


(3.6 ± 2.1) X 10"® 


190 


276 


216 


77 


14 


2 


g/cm^ 


(8.7 ± 1.4) X 10" 


* (2.69 ± 0.45) X 10"* 


(4.86 ± 0.93) X 10"® 


(6.3 ± 1.9) X 10"® 


(2.7 ± 1.9) X 10"® 


298 


151 


126 


40 


10 


2 


g/cm^ 


(9.9 ± 1.4) X 10" 


* (3.25 ± 0.48) X 10"* 


(5.2 ±1.0) X 10"® 


(8.8 ± 3.0) X 10"® 


(4.9 ± 3.5) X 10"® 


499 


272 


247 


95 


30 


4 


g/cm^ 


(6.61 ±0.87) X 10 


(2.35 ± 0.31) X 10"* 


(4.57 ±0.71) X 10"® 


(9.9 ±2.1) X 10~® 


(3.5 ±1.8) X 10"® 


957 


723 


735 


380 


124 


13 


g/cm^ 


(3.79 ± 0.45) X 10 


(1.49 ±0.17) X 10"'' 


(4.11 ± 0.50) X 10"® 


(8.6 ± 1.2) X 10"® 


(2.34 ± 0.70) X 10"® 


TABLE V. Near vertical neutrino intensities at 960 g/ 


'cm^ at zero geomagnetic cutoff, for ID and 3D TARGET calculations. 


All intensities are in (cm^ s i 


sr GeV/c)-^ 








Momentum 


(I'm +1^^) 


(i^M ± I'm) 


(i^e ± Ve) 


(z/e + Ve) 


(GeV/c) ID Intensity 


3D Intensity 


ID Intensity 


3D Intensity 


0.36 




0.32 


0.31 


0.15 


0.14 


0.71 




0.063 


0.062 


0.027 


0.026 


0.89 




0.035 


0.035 


0.014 


0.014 


1.4 




0.0099 


0.0099 


0.0038 


0.0038 


1.8 




0.0051 


0.0052 


0.0019 


0.0019 


2.8 




0.0014 


0.0014 


4.1 X 10"" 


4.3 X 10"" 


5.6 




1.8 X 10"'' 


1.8 X 10"" 


4.0 X 10"® 


4.2 X 10"® 
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2.3 X 10"® 


2.3 X 10"® 


3.7 X 10"® 


3.6 X 10"® 


22 




2.6 X 10"® 


2.6 X 10"® 


2.9 X lO"'' 


3.0 X 10"'' 
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